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The axial stiffness of MWCNTs is demonstrated to be determined only by several external shells 
(usually 3-5 and up to 15 for the extremely large nanotubes and high elongations) what is in a good 
agreement with experimentally observed inverse relation between the radius and Young modulus 
(i.e., stiffness) of MWCNTs. This result is a consequence of the van der Waals intershell interaction. 
The interpolating formula is obtained for the actual axial stiffness of MWCNT as a function of the 
tube external radius and elongation. 

PACS numbers: 61.48.De,62.25.-g,62.20.de 



A. Introduction 

The unusual and even sometimes wonderful mechanical properties of carbon nanotubes and their bundles make it 
possible to use them just now for the wide range of applications. As an example, nanotubes can act as an reinforcement 
of different materials (plastics, hydrocarbon resins, nanocomposites etc.) where their extremal bending flexibility and 
axial stiffness are of a great interest. [IH3] 

It is established that SWCNTs can sustain strains larger than 10% of tensile deformation prior to fracture [U [5] 
and their deformation is completely reversible (i.e., elastic) subjected to strains of more than 4%.[6H9] A lot of works 
were dedicated to investigation of elastic properties of single-walled nanotubes (SWCNT) and multiwalled nanotubes 
(MWCNT). [5"1I5HTE] Theoretically MWCNTs were studied for uniform axial stresses at both of their ends, [15] but 
this is not the only possible loading type. Apparently in most cases of axial tension only the external shell of MWCNT 
is affected by the imposed load and internal shells are involved into the considered process only due to van der Waals 
intershell interaction (the simplest case of such situation is the deformation of the capped MWCNT). 

It is obvious that in the case of uniformly deformed n-walled CNT the intershell distances Ar^ — ri~ri + \ = {l—ve)do 
remain equal (where i = 1,2, . . . ,n— 1 and i = 1 corresponds to the outermost shell), but when the load is imposed 
only on the external shell the distances between other shells are not equal and increase with i (for e > 0). Due 
to strong nonlinear dependence of the intershell interaction on the intershell distance the difference in deformation 
energy between both mentioned cases for certain specific elongation may be significant and should be studied. 

Some experiments showed that the effective Young modulus of the MWCNT is inversely proportional to its radius 
(for radii in range from 4 nm to 20 nm) [20] and such results contradict the assumption of the equal deformations 
of shells. Some theoretical works also showed the dependense of the effective Young modulus on the number of wall 
[TH1I2I] while other did not show any dependence on radius even for SWCNTs. [191 12"2"] 

In this paper we use the stiffness k ~ Yd instead of Young modulus Y in order to avoid the uncertain parameter 
"wall thickness" d. Some authors stand on using this parameter because when considering flexural deformations of 
nanotube (within the string approximation) we should work with Yd 3 but it is evident that we always can switch 
from one pair of independent parameters to another such pair (for example, from Y and d to k and 7, where 7 denotes 
some flexural characteristic). 

B. Axial Stiffness of SWCNT 

Since the investigation of stiffness (or Young modulus) of single-walled nanotubes is not the main objective of this 
work we will describe it briefly presenting here without derivation only most important for our purposes expressions. 
In the simplest case the axial stiffness of a cylindric shell (such as a nanotube) with surface density m^r (m is the 
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mass of each atom and r is their quantity per unit area) is defined by the expression 

R R 

k(R, L) = 2irm,QTC 2 — —a — , (1) 
L L 

where R and L are the shell radius and length, c is the longitudinal sound velocity and a is the constant parameter 
for graphene-based structures. Certainly previous expression is correct for nanotubcs of different radii only if shells 
properties are not radius-dependent. In the case of nanotubes with extremely small radii (R ~ do where do = 0.34 
nm is distance between graphitic planes conditioned by the van der Waals interaction) the re-hybridization of atomic 
orbitals lead to perceptible changes of its mechanical properties (the raising of effective Young modulus is observed 
on quantum dynamics simulations for SWCNTs with R < 6 nm |231 I24|). but for tubes with radii more than one 
nanometer such changes may be ignored. 

The sound velocity c may be obtained from the microscopic models of nanotube (or graphene) as a velocity of 
acoustic phonons. Further we will use the value c= 18.4 km/sec [35] which lead to aw 1632 kg sec -2 . If necessary 
the expressions for effective Young modulus and Poisson ratio also may be derived by comparing equations of motion 
from both microscopic and continuum models. 

By ([I]) one can calculate the idealized Young modulus of SWCNT 

Y = k -s=4 = a Ai~ o - 73TPa > 

where S is the effective surface of the nanotube cross section and the commonly accepted value of the "wall thickness" 
d = c?o is used. On the other hand the experimentally measurable effective Young modulus of the macroscopic bundle 
of SWCNTs with equal radii R depend not on the effective surface of its cross section but on the total cross section 
surface of the bundle Sb ~ Nn(R + do/2) 2 (where N is the number of tubes in the bundle): 



S b S b Tr(R + d a /2) 



a 1 

7T R' 



One can see that the effective Young modulus of such a bundle should be inversely proportional to the nanotubes 
average radius. 



C. Van der Waals intershell interaction 




FIG. 1: The definition of the main parameters of the considered DWCNT. 



Within the continuum approximation the van der Waals intershell interaction energy for the DWCNT depends only 
on shells radii (r 1; r 2 ) and lengthes (L 1 ,L 2 ) [Fig- 0- hi terms of hypergeometric functions it can be expressed as 
follows (assuming that L\ > L^) 

, 3 /21 <P(^,R 1 ,R 2 ) $ (%,R 1 ,R 2 )\ 
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where 76 = 2.43 x 10 24 J-nm 6 and 712 = 3.859 x 10 27 J ■ nm 12 are attractive and repulsive constants of the 
Lennard- Jones potential, [27J r is the surface density of carbon atoms and 

1 4i?ii? 2 \ + i? 2 ) 2,/ r d0 



$ (J, J? l5 Jfe) .- 2 Fi ' !' (j^ J ~ 2tt /_ 



+ i?| - 2R 1 R 2 cose) J 

In the case of unstrained graphene (or nanotube) r = r = — ^ fc2 =38.2 nm -2 (6 = 0.142 nm). 

The length and radius of the nanotube under the tension take the values l = L (l + e) and r = i? (l — fe) where e is 
the specific elongation of nanotube, v — 0.17 [2"8"H3Ti] is the graphene Poisson ratio, L and R are length and radius 
of unstrained nanotube. In this case the nanotube surface S and the surface density of atoms change correspondingly: 

Sj = 2TTRjLj(l - U£j)(l + £j), Tj - 



(1 + Ej )(1-i/ £j ) 

The interaction energy ([2| for the case of strained shells (ei and £2 respectively) takes on form 



(3) 



U{Ri,L 1 ,£i,R 2 ,L 2 ,e 2 ) = §7rViT 2 rir 2 Z 2 {^112 ( ri 2 + V 2 )n "76 (ri+r2 )5 

3„3 1 Jd d r ^21„ g (^^1(1-^1)^2(1-^2)) $(§ .^(l-^Q.jfaCl-vea)) ^ 

2 W l+ £l T ^1^2^2^32 712 (fl^i-^e^ + ^fi-^))!! 76 (i_„ ei )+fl, 2 ))5 ^ 



From now on we will denote by Uw the van der Waals deformation energy that is the contribution of the intershell 
interaction into the total deformation energy of double-walled nanotube 

AC/V(i?i, Li, ex, R 2 , L2, £2) := U(Ri,Li,ei,R 2 ,L 2 ,s 2 ) — Uo(Ri, L\, R 2 , L 2 ). (4) 

D. Axial Stiffness of DWCNT 

The total deformation energy of DWCNT consists from the deformation energies of both shells and the additional 
(to unstrained state) intershell interaction energy: 



E{e\, £2, Ri, R2, Li, L 2 ) = kl 2 Ll + fc ^L 2 -f AUw(Ri, Li, £1, i?2, L 2 , £2) — 

^§^- + ^§^ + AU w (R 1 ,L 1 ,e 1 ,R 2 ,L 2 ,e 2 ) = ( 5 ) 
§ (RiLief + R 2 L 2 £ 2 V ) + AUw{Ri, Li, £1, R 2 , L 2 , e 2 ). 



If the both shells are of equal length (L\ =L 2 =L), the total deformation energy of DWCNT will be linear in its 
length and the inner shell relative deformation will depend only on that of the outer shell and shells radii. 

We assume that the external force acts only on the outer shell and as a result the DWCNT axial stiffness should 
be introduced as 



L{ d£{ 



iB(..,es,..) 



Obviously, £2 < £1 when R\v £\ < do, and for the fixed value of £\ the £2 will increase strongly with increasing of 
DWCNT radius [Fig. |]. 

E. Axial Stiffness of MWCNT 

Like that of DWCNT the total deformation energy of TV-walled MWCNT is determined by the following expression 

N JV-1 

E(e,R,L,N) = ^J2 R i L i £ i + J2 AU w (Ri,Li,e i ,Ri + i,L i+1 ,e i+1 ), (6) 

»=1 i=l 

where e = (si, £jv)» R = (Ri, Rn) and L = (Li, Ljv). 
The axial stiffness in this case is defined as follows 

1 u R T m 1 d 2 E(e mm ,R,L,N) 



L\ 3e\ 



(7) 

E(e m i„,...)=min c _ ..., e B(e,R,L,JV) 
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Unlike the DWCNT all but innermost shells of MWCNT have the following neighbor which hampers the shell's 
deformation of comparing to the case of DWNCT so that the relative deformation of inner shells e\ decrease rapidly 
with i. As a result the actual axial stiffness of MWCNT grows with number of shells much slower than its "ideal" 
stiffness kideal — Tl,k{Ri,Li). Actually for MWCNTs with external radius R < 10 nm and arbitrarily number of 
shells under the strain with e\ < 0.05 only four external shells contribute to the total stiffness and this number grows 
slowly up to s» 10 for the nanotubes with R > 25 nm [Fig. [3]. 

For 7V-walled MWCNTs with external radii more than 3 nm (assuming that N is greater than the number of shells 
actually involved into the deformation) the stiffness can be fitted by the following expression with an accuracy of 1-3 
per cents (for < e± < 0.1) 



. 10~ 7 
k(R, L, e) = — - — exp 



2.3 



R-d, 



0.267 



1+e 



26.2-44 



R 



0.276 



(8) 



F. Discussion 



The analysis of axial stiffness of the ideal multiwalled nanotubes based on the van der Waals intershell interaction 
shows that only several MWCNT's external shells contribute to its total stiffness under the load imposed only on the 
external shell. The number of contributing shells in fact less than 5 for nanotubes with R < 10 nm and never exceed 
15 even for extremely thick nanotubes (figure [3]). As a result the MWCNT actual stiffness may be several times lower 
than its ideal stiffness which is obtained under the assumption that all shells are loaded evenly. This fact should be 
taken into account while MWCNT-based towlines, cables and armoring elements are constructed. 

Contrary to single-walled nanotubes, which are at least for small elongations subordinated to the Hooks law F — 
keL, k = const, the stiffness of MWCNTs even for rather weak loads depends linearly on the tube elongation (k ~ e) 
as a result of the load-induced gradual involvement of inner shells. 

All calculations presented here in support of our assertion should be performed for particular values of the Poisson 
ratio and sound velocity, which are (in the case of ideal lattice) uniquely determined only by the lattice structure 
and interatomic force constants. However up to now there is a large scatter of values of Poisson ratio obtained by 
different authors [y = 0.14 - 0.19, gS] v = 0.17 - 0.22, [3D] v = 0.19,0 |3T] v = 0.19 - 0.22, [3J v = 0.20 - 0.23, [33J 
v = 0.25 - 0.27, [33] v = 0.27 - 0.28 QJ]). In most cases v G 0.16 - 0.3. Here we used the value v = 0.17 which is a 
mean value of Poisson ratio for finite graphene sheets of different sizes [30] and chiral SWCNTs.[30] It is also close to 
the corresponding magnitude along the basal plane in graphite v = 0.16. [28] 
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FIG. 2: The dependence between the outer and inner shell relative deformations (ei and £2 correspondingly) for DWCNTs 
of different radii and unstrained intershell distance do. For very small deformations the behavior of a DWCNT is determined 
by its actual unstrained intershell distance which is not likely to be equal to the energetically optimal distance do due to the 
nanotubc's radius discrete nature. 
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In view of the existing data scattering we have looked how the calculated MWCNT stiffness depends on the used 
values of Poisson ratio and sound velocity. It clear that for Poisson ratios higher than v = 0.17 the strain transmission 
from the outermost shell into the depth of the tube raises what leads to higher values of nanotube stiffness. However 
even for v = 0.3 (which exceeds most of its theoretical and experimental estimates for SWCNTs) it appears that 
only 12 — 15 external shells contribute to the k even for MWCNTs with R > 25 nm and e > 0.05. In other words, 
despite the fact that actual stiffness of MWCNT grows almost linearly with increasing of the shells' Poisson ratio 
moving to its "ideal" value, nevertheless the latter remains many times higher. ( The values of axial stiffness for 
varying values of Poisson ratio may be found with an accuracy of several per cents using the interpolating formula 



The sound velocity actually depends on the quality of the shell's lattice and connected with the stiffness of each 
shell and "ideal" total stiffness of MWCNT by the relation k ~ c 2 . It appears that the sound velocity and the actual 
stiffness of MWCNT are linearly dependent in the interval c = 17000 — 23000 km/ sec as a consequence of a weakening 
of inner shells strains due to the preferred accumulation of the van der Waals interaction energy in subsurface shells. 
So if used here c= 18.4 km/sec is underestimated the number of shells actually contributing to the MWCNT stiffness 
would be even lower than indicated above. 

Note that the stiffness of perfect MWCNTs is perceptibly greater than that of MWCNTs with rather low defect 
concentration, [35 but their ratios to the corresponding ideal values are in the opposite relation (i.e., the actual and 
ideal rigidities of some MWCNT would be most close in values in the case of pretty defect MWCNT). 

It also should be noted, that some of MWCNT inner shells may be segmented. In accordance to our analysis such 
situations do not affect substantially the axial stiffness in all cases when the cumulative length of "broken" shell is 
close to length of its neighbors and gaps between shell's parts are about the range of van der Waals interaction (~ 1 
nm), that is the elastic properties of MWCNT are resistant to fracture and minor damage of some of its inner shells. 

In that case when one of inner shells with number i appears to be sufficiently shorter than neighboring shells the 
MWCNT can be considered as consistent of two placed in series parts with different number of walls and as a result 
different rigidities (the first part of length L\ — Li is a (i — l)-walled MWCNT and the second part of length Li is a 
n- walled MWCNT). The effective stiffness for such a series is k = , where k\ and k^ are determined by &7s. 

Thus for applications where the highest possible axial stiffness of single nanotube or nanotube bundle is needed the 
4-5- walled nanotubes of minimal diameter are quite sufficient. 



k{v) « fcM i + 
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FIG. 3: The dependence of MWCNT stiffness on the tube's external radius R\ and the number of its shells (in brackets) for 
the £i = 0.05 and L — 100 nm. The asymptotic value of the number of shells involved into the deformation is about 15 even 
for the extremely large MWCNTs (for tubes with R < 10 nm there only four shells are actually involved). 
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